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1. INTRODUCTION 

Each distribution system has many literals, load taps, balanced and non-balanced loads, and 
different types of conductors, which complicate the calculation of fault location in this system. Many studies 
have been done to locate faults in distribution systems [1]-[5]. In [1], short line model and phase components 
were used for locating. The sensitivity of this method to the fault resistance is high. In [2]-[4], the distributed 
parameter line model and the main component of voltage and current are used. In These articles the presence 
of literals has not been considered. In [5], the distributed parameter line model has been used for a single- 
phase-to-ground fault, but the existence of literals and unbalanced performance of the network have not been 
analysed. 

Permanent high impedance faults in distribution networks are generally due to their low level of 
voltage. The low voltage of the distribution network and high ground resistance cause low fault current in 
this type of faults. Therefore, protections in the network, such as distance or overcurrent relays, will not be 
able to identify them [6]. These types of faults generally occur in two ways, such that the conductor contacts 
a high impedance object, such as branches of the tree, or the conductor is cut off on the ground [7]. To date, 
most of the work in the distribution system has been limited to detecting high impedance faults, these include 
the use of numerical fault analysis [8], neural network [9] and the presentation of a wavelet transform 
algorithm [10]. Regarding the location of high impedance fault in the distribution system, one can refer to the 
use of wavelet transform in [11]. Since it is impossible to ignore the capacitive effect of the line during the 
fault in the distribution system, modeling this effect is very important for having high accuracy in obtaining 
the fault location. In recent years with the development of telecommunication systems, phasor measurement 
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units (PMU) have changed to one of the most important systems for Wide Area Monitoring, Protection and 
Control (WAMPC) of power systems [12]-[14]. These equipment, using strong signal processing technology 
and a global positioning system (GPS), enables to simultaneously calculate the voltage and current of the 
branches on which the PMU is installed [15]. In [16] and [17], PMU placement in the distribution network is 
proposed to estimate the network state, and furthermore, this technology is used in fault locating of power 
systems [18]-[20]. 

In this paper, considering the importance of the capacitive effect of line, an algorithm is proposed to 
improve the accuracy of impedance-based fault locating methods, using the distributed parameter line model 
in phase domain, which, contrary to previous methods [1]-[6], has lower sensitivity to resistance, location 
and type of the fault. In the proposed method, using phasor measurement units installed in the distribution 
system, voltage and current of all parts of the distribution system is calculated. Then, two buses where the 
fault occurred in the line between them is determined and with the other calculations, the fault location is 
computed. The proposed method is evaluated on 25 kV network (Saskpower) in Canada [21] in the presence 
of literals and static loads. The results showed high accuracy of the algorithm. 

In the remaining of the paper, the second part describes the proposed method for calculating the 
fault location and determining the type of fault in the line in which the error occurred. The third section 
describes the PMUs placement and how to apply the proposed method to the distribution networks. In the 
fourth section, the simulations and the results for the evaluation of the proposed method have been 
investigated and, finally, a conclusion has been made. 


2. PROPOSED METHOD 
2.1. Fault classification 

In this section, considering the distributed parameter line model in the phase-domain and using 
impedance-based method, the article sets to calculate the fault location, type (phase-to-phase and phase-to- 
ground) and resistance. One section with Ax length of the line can be shown in Figure 1. According to Figure 
1, the Equation 1 can be written as follows: 


v 
; (1) 


where: AV: Voltage drop on series impedance of elemental length of the line 
AI: Capacitance current of elemental length of the line 
z: series impedance matrix of the line per km 
y: shunt admittance matrix of the line per km 


— 
I I+AI 
> —> 
L.Ax R.Ax | 


Figure 1. Elemental section of line with length of Ax 


The differentiation of Equation 1 with respect to x leads to Equation 2. In power lines, phases 
mutually affect each other due to mutual coupling. Using the method of symmetrical components (matrix 
(3)), independent phase components and Equations 4 and 5 can be obtained: 

2 

z= ZV 


H = yzl 


(2) 
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Solving Equations 4 and 5 and using Equation 6, the sequence of the voltage and current of each 
section can be calculated. K represents the sequences of positive, negative and zero. 
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where: l: the total length of section 
VR: voltage of end bus of section 
IR: current that is inputted to end bus of section 
As will show in subsection 3.2, using phasor measurement units installed in the distribution system, 
the voltage and current of the two S and R buses, the fault occurred between them, can be calculated. In 
Figure 2 Using Equations 7 and 8, the voltage and current on the two sides of the fault point can be obtained. 
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Figure 2. Line between S and R buses that fault is occurred 
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where: VŠ : Voltage of fault point from sending side 
Vee : Voltage of fault point from receiving side 
I 5 : Current of fault point from sending side 


I z : Current of fault point from receiving side 
Since the voltage is equal at the point of fault, the following equations can be written: 


veya (9) 
Vor = Vep (10) 
Vor = Vep (11) 


The only unknown factor in the above equations is the fault location per line length unit (D) that can 
be obtained through solving Equations 9, 10 and 11 using Newton-Raphson algorithm. When the network is 
in the form of single-phase lines, instead of Equation 6, the following equation is used: 


Gcosh(g 1) z, sinh(g 1) Ù, 


Ê. 
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In the next sub-section, the fault type detection method in the two groups of phase-to-phase and 
phase-to-ground is separately determined and the fault resistance will be calculated. 


2.2. Fault Type Determination and Fault Resistance Calculation 
To determining the type of fault, using transfer matrix and according to the equations 13 and 14 the 


phase voltage and current at the point of fault can be estimated: 


yu é 1 1 tee 
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2.2.1. Phase to Phase Fault 
For example, the phase to phase fault for A and B phases by impedance fault ZFab is shown in 


Figure 3. The fault current can be obtained as follows: 


Lpa = eer E l receive (15) 
Tp = I, send ~ I, receive (16) 
Tn = Lesena E L receive (17) 


Int J Elec & Comp Eng, Vol. 8, No. 5, October 2018 : 2709 - 2720 


Int J Elec & Comp Eng ISSN: 2088-8708 O 2713 


If one of the equations 18, 19 and 20 holds true, there is phase to phase fault: 


Tig any Tn (18) 
Lia => Le (19) 
L =e Tne (20) 


Va- V, 
Z = Fa Fb (2 1) 


Figure 3. Phase to phase fault model 


2.2.2. Phase to Ground Fault 
If Equations 18, 19 and 20 do not confirm, it can be concluded that a phase to ground fault has been 


occurred. This fault is modelled as shown in Figure 4. The fault impedance matrix can be written like 
equation 22: 


yy Ù A + oe, Lig h Ug U 
e"u & L 
a ay AtA a (22) 
eru T L 
MA Z, Zn. + Z Bark 


As a result, due to the fact that the phase voltage and current of the fault point are known, the 
impedance matrix can be calculated and according to the values obtained for the fault impedance of equation 
22 and the fault current, one can determine which types of faults (Single-phase to ground, two-phase to 
ground and three-phase to ground) occurred in the distribution system. 


S Zfg 


Figure 4. Phase to ground fault model 
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3. IMPLEMENTATION OF THE PROPOSED METHOD ON DISTRIBUTION SYSTEMS 

The proposed method for fault location is based on assumption that the voltage and current of the 
two buses in which the fault occurred in the line between them are known. However, only the voltage and 
current of a limited number of buses on which PMUs are installed are available. Therefore, in the rest of the 
article, a method for PMUs placement with respect to the proposed algorithm is presented to obtain the data 
of other buses through using the PMUs data. 


3.1. PMU Allocation in Distribution Systems 
3.1.1. Network Observability Using PMU 

PMU installed on a specific bus can measure the voltage magnitude and phase angle of bus as well 
as the magnitude and phase angle of branches currents attached to that bus. As a result, the magnitude and 
phase angle of voltage and current of buses connected to the bus that equipped with PMU can also be 
estimated using basic rules (KVL, KCL) of electricity. Consider, for example, Figure 5, the buses in which 
the PMU is installed (PMU buses) are direct observable, and those buses that are connected to the PMU 
buses are indirect observable. The others are not observable [22-23]. 


Not Observable network e- ~ Observable network — 


Not Observable pat —— — = 


Figure 5. Network observability using PMU 


3.1.2. PMU Allocation Algorithm 
PMU placement should provide the following conditions: 

a. According to the placement of PMUs the voltage and current of all buses from two sides can be 
calculated. This means that when fault is occurred in one side of a bus, the data of bus can be acquired 
from the other side that is intact (completely observability). 

b. In all buses not connected from one side to the other bus (end buses), a PMU is installed so that 
information is available in the event of a fault on the other side of the end bus. 

After allocating PMUs with considering above conditions, the voltage and current of all buses can be 


calculated at any time. 


3.2. Feeder’s Data Calculation 

To better explain the contents of this subdivision, for example, a single-line diagram of a sample 
distribution system is shown in Figure 6. According to PMU placement algorithm the PMUs are installed on 
buses 1, 5, 6. Therefore, the voltage and current of these buses will be simultaneously and directly available 
for fault locator. The calculation of the voltage and current of buses 2 and 4 in the normal mode of the 
network will be as follows: 

It is assumed that loads in the buses where PMU is not installed are static. From bus 6, using 


Equation 6, voltage and current I 
23: 


Væ Can be obtained. Therefore I, 


aae Nate can be calculated using Equation 


be 
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V and from bus 5, I 


ma V V,,. are estimated according to equation 6 and then 


35abc? 3abc 


is obtained via (24): 


V, 
I = I + 3abc (24) 


32 abe 35 abe 


s2abc > the current I,,,,. is estimated, and: 


_ I = I = V, abc (2 5) 
24abe ~~ 21abc 23abe Z 
L2 
Thus, Viste 2 Vo abe ? Lave ? Dysabe are known. 


Figure 6. A single-line diagram of a sample distribution system 


Similarly, according to the proposed method for installing PMUs in each distribution system, the 


calculation of the voltage and current of buses with the algorithm presented in this subsection is 
generalizable. In the event of a fault in the network, voltages and currents are also calculated, except that 
symmetric components are used for observation. 


3.3. Fault Location Algorithm 


a. 
b. 


The proposed fault location algorithm is presented as follows 
The synchronous data are received from all PMUs in a cycle. 
The calculation of the voltage and current of buses of the entire distribution system according to the 
algorithm provided in section 3.2 (estimation of voltage and current of all the pairs of buses at any time) 
Equation 6 should be checked for all pairs of buses. Pair of buses in which Equation 6 does not hold are 
considered as sections with fault. (Go to the next step). If Equation 6 holds for all pairs of buses, go to 
step1. 
According to Equations 7 to 11 and using Newton-Raphson's algorithm, fault location in line unit length 
(D) is estimated for all pairs of buses that are considered as the section with fault. 
Fault has occurred between two buses where 0 <D <1 and D is the fault location per line length unit. 


V_,. are estimated according to the equations (13) to (17). The 


The voltage and current of fault point I, V,.,,. 


type of fault and fault resistance can be determined according to sections 2.2.1 and 2.2.2. The algorithm 
flowchart is shown in Figure 7. 
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Receive Data from PMU 


Calculate current of fault point 
(one cycle) 


according to equation (7) & (8) <———_ 
Isp,Ipr 


Calculate voltage and current of all 
the pair nodes in Distribution network 
VS,VR,IS,IR 


elect two nodes Ifa=Iasp-Iapr 


(respectively) Ifb=Ibsp-Ibpr 
Ifc=Icsp-Icpr 
Equation 6 is hold for two i 
YES Re 3 
nodes? is investigate? Ifa=-Ifb?? 
| Or 
NO Ifb=-Ifc?? 
t - Or 
Calculate voltage of fault point NO Ifa=-Ifc?? YES 
according equation (7) & (8) Phase to Phase to 
Vsp,Vpr ground Fault 4- —* Phase Fault 
is occured is occured 
Solve the equation Vsp-Vpr=0 using } 
Newton Rawson algorithm and obtain} NO 
unknown D The fault resistance is obtain The fault resistance is 
from Eq. 22 obtain from Eq. 21 
YES 


YES 
D is location fault per 
length 


Figure 7. Flowchart of proposed method 


4. SIMULATION RESULTS AND ANALYSIS 

To evaluate the proposed method using a real system, 25 kV distribution system of Saskpower in 
Canada [21] was studied. The 25 kV single-line diagram is shown in Figure 8. The implementation of the 
proposed algorithms was done in MATLAB and then the simulation of 25KV distribution system was carry 
out by EMTP-RV software. In the conducted simulations, distributed line parameter model was used. All 
loads were considered as static load. The percentage of error in fault locating is obtained from the following 
equation: 


actual fault location-estimate fault location , 
percentage of error= —————————— l 100 (26) 


length of line 
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Figure 8. Diagram of 25 kV power distribution system of Canada 


Various simulations for three different pairs of buses with different types of faults (single-phase to 
ground, phase to phase, 3-phase to ground) were done on the 25 kV distribution system. The results of the 
implementation of the proposed method are shown in Table 1. It should be noted that the line between buses 
(3-4) and (6-13) is three-phase, and the line between buses (18-20) is single-phase. 

To demonstrate the effectiveness of the proposed method for faults with low fault impedances, 
faults with 10 ohms of resistance were simulated, and the proposed method in these cases were evaluated, 
and the results are shown in Table 1. In this table, in addition to the fault location estimation, the fault 
resistance was also estimated. The results show the high accuracy of the proposed method. 

As shown in Table 1, fault locating is implemented between three pairs of buses in Canada 25 KV 
distribution system. Faults have occurred at three distinct distances between two buses. In order to 
demonstrate the efficiency of the proposed method in various fault resistances, faults have been investigated 
in three distinct fault resistances (10, 100 and 300 ohms). The maximum error rate in the fault location is less 
than 0.15%, which indicates the high accuracy of the proposed method. Also, regard to the numbers obtained 
in different conditions, it can be concluded that the accuracy of the method does not depend on the fault 
conditions and the results indicate the low sensitivity of the proposed method toward the fault resistance, 
fault location, and fault inception. Figure 9 to 10 shows fault percentage of error in relation to the fault 
resistance of 10, 100 and 300 ohms. As shown in figures, fault locating for each position in different fault 
resistance is close to each other. 
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Figure 9. (a) Percentage of error between bus 3 and 4 considering different conditions. (b) Percentage of error 
between bus 6 and 13 considering different condition 
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Figure 10. Percentage of error between bus 18 and 20 considering different conditions 


Table 1. Simulation Results 


Fault Actual fault Type of fault 
resistance location Single phase to ground Phase to phase Three phase to ground 
Fault Fault Percentage of Fault Percentage of Fault Percentage of Fault 
between Bus Location Error Resistance Error Resistance(ohm) Error Resistance 
(3-4) from Bus 3 (location) (ohm) (location) (location) (ohm) 
10 ohm 1005.75 0.0082190 9.9861 0.01094 10.0134 0.01456 9.9921 
100 ohm meter 0.0074577 100.0379 0.007523 100.0685 0.007612 100.052 
300 ohm 0.0074354 300.178 0.007487 300.2094 0.007462 300.1934 
10 ohm 2212.65 0.0070608 9.9868 0.01047 10.0069 0.015064 9.9924 
100 ohm aer 0.065043 100.016 0.006545 100.0382 0.006605 100.0279 
300 ohm 0.006525 300.1001 0.006526 300.1224 0.00652 300.1126 
10 ohm 2816.1 0.006896 9.9896 0.01363 10.0042 0.02107 9.994 
100 ohm pei 0.0029392 100.0084 0.0032659 100.0245 0.003713 100.0175 
300 ohm 0.0028556 300.0644 0.002908 300.0806 0.002971 300.0739 
10 ohm 603.5 0.015352 10 0.031029 10 0.047291 10 
100 ohm sate 0.0031594 100.0002 0.0043964 100.0002 0.0059012 100.001 
300 ohm 0.0026904 300.0006 0.0028942 300.0006 0.0031919 300.003 
10 ohm 1327.7 0.014013 10 0.028465 10.0001 0.043621 9.9999 
100 ohm neter 0.0021049 100.0001 0.0035189 100.0006 0.0050663 100.0001 
300 ohm 0.0014482 300.0003 0.0017519 300.0004 0.0021499 300.0004 
10 ohm 1689.8 0.012104 10 0.024591 10.0001 0.037811 10 
100 ohm méter 0.0015841 100.0001 0.0029080 100.0004 0.0042982 100.0004 
300 ohm 0.00087521 300.0002 0.0012204 300.0012 0.0016305 300.0012 
10 ohm 804.75 0.080157 10.1206 - - - - 
100 ohm PAo 0.028213 107.1013 
300 ohm 0.023404 298.0522 - - - - 
10 ohm 0.038865 9.998 
100 ohm P 0.022889 99.7842 - - - - 
300 ohm 0.01899 298.1764 
10 ohm 0.14 9.9981 - - - - 
100 ohm pecs 0.035857 99.7856 
300 ohm 0.019663 298.1952 - - - - 


In order to validate the performance of the proposed method, Table 2 examines the error rate in the 
fault locating in distribution systems with different methods. As shown in table 2, the method presented in 
this paper has a relative superiority to other methods. 


Table 2. Comparison between Different Methods 
Mean of error (%) 


Ref [24] Ref [25] Ref [26] 
Fault type Propose ameo Fault res. (10ohm) Fault res. (10 ohm) Fault res. (1000hm) 
Single phase to ground 0.021 - 2.293 1 
Phase to phase 0.009 0.065 - - 
Three phase 0.013 0.085 2.294 - 
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5. CONCLUSION 

Finding a fast and accurate fault locating method has been of particular importance. In the 
meantime, it should be considered that the location of the fault in the distribution system is more complex 
than the transmission system. Due to the increasing use of PMUs in power systems for network monitoring, 
in this article, the PMUs have been used to provide a method for fault locating. In the proposed method, first, 
according to the system topology, the location of the installation of the PMUs should be located so that the 
data of all buses can be calculated from both sides during the faults. If a fault occurs in the system, the 
information obtained from both sides of a bus is not consistent with each other. Subsequently, based on the 
proposed algorithm the fault location and type will be estimated. In the proposed method, according to the 
characteristics of the PMUs, the asymmetry and the presence of several laterals in lines, in the fault locating 
are ineffective. The maximum percentage of fault locating error is less than 0.15%, indicating the high 
accuracy of the proposed method. The results also showed the low sensitivity of the proposed method to 
resistance, location and type of the faults. Also, the evaluation of the proposed method and comparison with 
other previous work on fault locating in distribution systems indicates the superiority of the proposed 
method. 
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